#!/usr/bin/env python

'''
program plots phase diagram for ferromagnetic triangular transverse
'''

import matplotlib.pyplot

B = [0,    2,     3,    3.5,  4,    4.5,     4.72]
T = [3.62, 3.355, 2.87, 2.51, 1.95, 1.335,   0]

matplotlib.pyplot.xlabel("$B/J$")
matplotlib.pyplot.ylabel("$T/J$")
matplotlib.pyplot.plot(B, T, "*-")
matplotlib.pyplot.title('Phase diagram')
matplotlib.pyplot.show()